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ABSTRACT 

We combine data from the extremely deep Hubble Space Telescope U (F300W) image obtained using 
WFPC2 as part of the parallel observations of the Hubble Ultra Deep Field campaign, with BVi images from 
the Great Observatories Origins Deep Survey (GOODS) to identify a sample of Lyman break galaxies in the 
redshift range 2.0 < z < 3.5. We use recent stellar population synthesis models with a wide variety of ages, 
metallicities, redshifts, and dust content, and a detailed representation of the H i cosmic opacity as a function 
of redshift to model the colors of galaxies in our combination of WFPC2/ACS filters. Using these models, we 
derive improved color selection criteria that provide a clean selection of relatively unobscured, star forming 
galaxies in this redshift range. Our WFPC2/F300W image is the deepest image ever obtained at that wave- 
length. The 10(7 limiting magnitude measured over 0.2 arcsec 2 is 27.5 magnitudes in the WFPC2/F300W 
image, about 0.5 magnitudes deeper than the F300W image in the Hubble Deep Field (HDF)-N. This extra 
depth relative to the HDFs allows us to directly probe the luminosity function about 0.5 magnitudes deeper 
than the depth accessible with the HDF data along an independent line of sight. Our sample of star-forming 
galaxies with 2.0 < z < 3.5 includes 125 objects, the majority of which show clumpy morphologies. We mea- 
sure a star formation rate density of 0.18 M0yr _1 Mpc , marginally higher than the value measured for the 
Hubble Deep Fields. 

Subject headings: cosmology: observations — galaxies: distances and redshifts 
galaxies: statistics — ultraviolet: galaxies 
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1. INTRODUCTION 

Great effort has been expended in the last decade to identify 
galaxies at high redshift. It is hoped that measurements of 
their luminosity and color evolution will place constraints on 
the star formation history of the universe and indirectly on 
the formation of large scale structure. In this search for high 
redshift galaxies, the Lyman continuum break technique has 
been successfully applied to detect galaxies at cosmological 
distances. 

A characteristic feature in the UV spectrum of star form- 
ing galaxies is the Lyman break, the sharp continuum discon- 
tinuity at 912 A. The break is caused by absorption of UV 
photons with wavelength A < 912 A in the outer atmospheres 
of massive stars that produce such photons, in the interstel- 
lar H i gas, and by intervening H i gas along the line of sight 
from the star forming galaxy to the observer. At a redshift of 
z ~ 3, the Lyman break is redshifted longward of the U filter, 
so that galaxies have very low (or undetectable) flux in U. If 
the galaxy is starforming, it will have a high flux redward of 
the Lyman break, causing it to be bright in observed B and 
blue in colors such as B - V or V - 1. A color selection that 
identifies "U dropouts", i.e., galaxies that are blue in B — V 
and very faint in U, would then preferentially select actively 
starforming galaxies in the redshift range 2 < z < 3.5. The 
same principle can be used to identify galaxies at higher red- 
shifts by moving to a redder set of filters. For example, the 
B dropouts which have z ~ 4.5 can be effectively identified 
using BVI data. 

Several research groups have identified and studied U 
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dropout galaxies (see Giavalisco (2002) for a review). Ground 
based studies, mostly carried out with the Keck telescope 
(e.g., Steidel et al. 2003 and references therein), have stud- 
ied samples with deep imaging and spectroscopy over rela- 
tively large areas of sky (typically about 1000 square arcmin), 
to a typical depth of Rab = 25.5 largely set by the spectro- 
scopic follow-up. These studies were complemented by sam- 
ples identified using the Hubble Deep Fields, two deep imag- 
ing surveys with the WFPC2 instrument aboard the Hubble 
Space Telescope (HST) (Williams et al. 1996; Casertano et al. 
2000, hereafter COO). The Hubble Deep Fields were used to 
identify and study U dropout galaxies (Madau etal. 1996; 
Dickinson (1998); COO) to much deeper limiting magnitudes 
(Rab ~ 28.0), but over a solid angle about two orders of mag- 
nitude smaller than the ground-based observations. The HDF 
samples were too faint for complete spectroscopic follow-up. 
The precision of HST photometry also ensures small random 
errors on color measurements. Another important advantage 
of the HST is its unprecedented angular resolution that allows 
for a study of the main morphological characteristics of the U 
dropout galaxies, even at high redshift. However, the small 
solid angle that they probe makes them susceptible to sta- 
tistical fluctuations caused by small number statistics, galaxy 
clustering and field-to-field variations. The amplitude of these 
fluctuations can be estimated by observing multiple fields sep- 
arated by a large angular distance. Presently, the two Hubble 
Deep Fields provide the only two such fields on the sky. 

The scientific success of the WFPC2 Hubble Deep fields 
motivated an even deeper set of observations with the Wide 
Field Channel of the Advanced Camera for Surveys (ACS) 
installed on the HST in March 2002. The Hubble Ultra Deep 
Field (HUDF; Beckwith etal. 2006) covered an area about 
twice that of each of the original deep fields and reaches about 
1 .5 magnitudes deeper in the i band. 

The ACS Wide Field Channel has been optimized for visi- 
ble and near-IR throughput and has very poor UV sensitivity. 



2 



Wadadekar et. al 



For this reason, the HUDF ACS observations were restricted 
to the BViz (F435W, F606W, F775W and F850LP) filters; U 
band data were not obtained. The HUDF thus provides the 
deepest available sample of B, V, and i dropout galaxies, but a 
study of U dropouts is not possible. However, in parallel with 
the HUDF ACS and NICMOS prime observations, data were 
gathered with all the other HST instruments in parallel mode. 
As part of these parallel observations, the WFPC2 instrument 
observed two fields using the F300W filter. Both these fields 
are considerably deeper than the HDF F300W observations. 
One of them overlaps with the coverage of the ACS GOODS- 
South field (Giavalisco et al. 2004) and thus has moderately 
deep BViz data from that survey. In this paper, we combine 
data from the GOODS survey with the WFPC2 parallel ob- 
servations of the HUDF to identify a new sample of U band 
dropout galaxies. Our deeper i/-band observations, relative 
to the Hubble Deep Fields, allow us to directly explore a re- 
gion of the luminosity function of U dropout galaxies that has 
never been directly probed before. 

Throughout this paper, unless otherwise stated, we use the 
standard concordance cosmology with VLm = 0.3, J7a = 0.7 and 
ftioo = 0.7. For magnitudes, we use the AB system of Oke and 
Gunn (1983). For convenience, we use the designation U, B, 
V and i while referring to the WFPC2/F300W, ACS/F435W, 
ACS/F606W and ACS/F775W filters respectively. 

This paper is organized as follows: in Section 2 we describe 
our data and data processing procedures. In Section 3 we de- 
scribe the simulations we carried out to identify the appropri- 
ate color selection criteria for our dropout candidates and how 
we applied the criteria to identify our dropout sample. Section 
4 describes the dropout morphology and Section 5 contains a 
discussion of dropout statistics. 

2. THE DATA 

In this paper we analyze the HUDF/WFPC2 (F300W) par- 
allels which overlaps substantially over the GOODS-S area 
i.e. the Epoch 1 -Orient 310/314 image (see Fig. 1 of de Mello 
et al. (2006) for a finding chart of this area relative to the 
GOODS-S field). The other WFPC2 parallel field (Epoch 2- 
Orient 40/44) overlaps with the Galaxy Evolution From Mor- 
phology and SEDs (GEMS; Rix etal. 2004) survey area; be- 
cause of the lack of B-band data, this region cannot be studied 
in the same way, and is not considered further in this paper. 
Each field includes several hundred exposures with a total ex- 
posure time of 323.1 ks and 278.9 ks respectively. Because 
of small changes in the position angle of HST during each 
observing campaign, data at each epoch cover 6.52 square ar- 
cmin, with 4.02 square arcmin at the full nominal depth (the 
three wide field chips from all datasets overlapping), and the 
remaining 2.50 between 25% and 75% of the observing time. 
About 87% (5.67 square arcmin) of the F300W image over- 
laps with the GOODS-S footprint, including almost the entire 
full depth area (see Figure^. 

2.1. Data Processing 

The data processing was carried out with the techniques 
developed for the WFPC2 Archival Parallels Project (APPP; 
Wadadekar etal. 2006), based on the drizzle approach 
(Fruchter & Hook 2002). We constructed a drizzled image 
with a pixel scale of 0.06 arcsec/pixel. The various proce- 
dures required to obtain a final combined image have been 
described in de Mello et al. (2006). For completeness, we 
repeat the data processing steps here. 

For Epoch 1 , a total of 409 WFPC2/F300W parallel images, 



with exposure times ranging from 700 seconds to 900 seconds 
were obtained in parallel with the prime ACS observations. 
Each of the datasets was obtained at one of two orientations 
of the telescope: (i) 304 images were obtained at Orient 314 
and (ii) 105 images were obtained at Orient 310. 

We downloaded all 409 datasets from the HST data archive 
along with the corresponding data quality files and flat field 
images. These datasets were all downloaded in a single re- 
quest to the archive to ensure that the "on the fly reprocessing" 
(OTFR) procedure was identical for all the images. World Co- 
ordinate System (WCS) coordinates listed in the FITS header 
of WFPC2 images may be inaccurate by as much as an arc- 
second, depending on the details of the HST guide star ac- 
quisition process. This inaccuracy can lead to relative offsets 
between images that contribute to the drizzle stack, as well as 
errors in the absolute positioning of the image on the sky. 

The processing pipeline of the APPP registers images rela- 
tive to each other using an automated procedure that employs 
matched source lists to compute the offsets. We modified this 
procedure by measuring the centroid positions of 4 predeter- 
mined stars in every image with respect to an arbitrarily cho- 
sen reference image. These stars are the only sources that are 
detectable in each individual exposure in our drizzle stack. 
The nominal WCS of each individual image in the drizzle 
stack was updated to reflect the improved relative alignment 
with respect to the reference image; after this WCS update the 
images were well registered with respect to the reference im- 
age. Given that we only had 4 stars to fit to, we were unable 
to determine independent rotation corrections to a sufficiently 
high precision, and therefore left unchanged the orientation of 
each individual image, correcting only for shifts in the posi- 
tion of the stars. 

Our procedure also accounted for the shifts in the relative 
positions of the WFPC2 detector as a function of time using 
the model of Casertano & Wiggs (2001). Geometric distortion 
coefficients determined by Kozhurina-Platais et al. (2003) for 
the F300W filter were incorporated into the procedure. Cos- 
mic ray rejection was achieved using the procedure outlined 
by Fruchter & Hook (2002). 

The processing was carried out in two stages. In the first 
stage, images with near 100% overlap were grouped together 
in groups of 20 images each. Each group was processed sep- 
arately. Such separate processing in groups were necessitated 
by memory limitations on our computer. The drizzled image 
stack of 20 images was used to identify cosmic rays and store 
that information in a cosmic ray flag image for each of the 20 
images combined. This step was repeated for all the groups. 
In the second stage, all 409 images were drizzled through to a 
final image using the cosmic ray flag images to reject cosmic 
rays. 

Offsets between our final drizzled image and the GOODS 
images were measured by matching sources in our image with 
the corresponding sources in the GOODS data, which were 
binned from their original scale of 0.03 arcsec/pixel to 0.06 
arcsec/pixel. Once the offsets between the drizzled WFPC2 
image and the GOODS images had been measured, all 409 
images were drizzled through again taking the offsets into ac- 
count, so that the final image was accurately aligned with the 
binned GOODS images. This final image had an identical 
World Coordinate System (WCS) as the binned GOODS im- 
ages; this means that all sources have identical pixel coordi- 
nates in all filters and can thus be used for accurate dual mode 
photometry, where the image in one filter is first used for de- 
tection and subsequently photometry is carried out in all filters 
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using the segmentation map of the detection image. The resid- 
ual scatter in source registration between the GOODS images 
and the WFPC2 image is random (20 milliarcsec r.m.s.) and 
is likely caused by the mismatch in the centroids of sources 
between the U and redder bands, which is probably due to the 
fact that the regions of peak UV emission in galaxies probed 
by the U band image may not coincide with the peaks of light 
emission at redder wavelengths, leading to random offsets in 
centroid positions. 

We weighted each drizzled dataset by its inverse variance 
map during the drizzle process. This map was computed ac- 
cording to the prescription of COO. This computation of the 
variance takes into account contributions to the noise from 
the sky background (as modulated by the fiat field), dark cur- 
rent and read noise. Contributions to shot noise from sources 
are not included. Each weight map is combined with masks 
that exclude (i.e., set to zero weight) pixels that are flagged 
as bad in the data quality files, and is then used to weight 
the combination of images in the drizzle process. The drizzle 
process also produces a weight map as one of its outputs. The 
weight map is an accurate measure of the effective exposure 
(and therefore depth) that is reached at any particular location 
in the image after correctly accounting for missing data due 
to cosmic rays, hot pixels etc. in each image that is being 
combined. 

The WFPC2 CCDs have a small but significant charge 
transfer efficiency problem (CTE) which causes some signal 
to be lost when charge is transferred down the chip during 
readout (Heyer, Biretta, et al. 2004). The extent of the CTE 
problem is a function of target counts, background light and 
epoch. Low background images (such as those in the F300W 
filter) at recent epochs are more severely affected. Not only 
sources, but also cosmic rays leave a significant CTE trail. 
We attempted to flag the CTE trails left by cosmic rays in the 
following manner: if a pixel was flagged as a cosmic ray, adja- 
cent pixels in the direction of readout (along the Y-axis of the 
chip) were also flagged as cosmic-ray affected. The number 
of pixels flagged depends on the position of the cosmic ray 
on the CCD, with more pixels flagged for higher row num- 
bers. Only one pixel was flagged for rows 1 through 100, 2 
pixels were flagged for rows 101 through 200 and so on. With 
this approach, we were able to eliminate most of the artifacts 
caused by cosmic rays in the final drizzled image. 

The lOcr limiting magnitude measured over 0.2 arcsec 2 is 
27.5 magnitudes over the ~4 square arcmin field that has 
all 409 datasets overlapping, which is about 0.5 magnitudes 
deeper than the F300W image in the HDF-N and 0.7 magni- 
tudes deeper than that in the HDF-S. Somewhat lower depth 
is reached in the regions where not all datasets overlap. 

3. MODELING THE COLORS OF HIGH REDSHIFT GALAXIES 

3.1. Evolutionary population synthesis models 

In order to identify galaxies in the redshift range of interest, 
candidates first need to be identified in an appropriate color- 
color space of the available data. Such a "selection function" 
may be affected by many factors: the particular color criteria 
adopted, the intrinsic dispersion in the UV spectral properties 
of the star forming galaxies (as influenced by differences in 
metallicity, stellar initial mass function, star formation history 
and dust properties), the cosmic variance along different lines 
of sight and the photometric measurement errors. This selec- 
tion function can be estimated from models or can be mea- 
sured directly, by determining enough spectroscopic redshifts 
to measure it empirically. Empirical determination for space 



based dropout samples is precluded for at least two reasons. 1 . 
faintness of the candidates makes spectroscopy difficult even 
with 10 m class telescopes and 2. clustering of Lyman break 
galaxies necessitates measurements over many independent 
sightlines to average over the effects of large scale structure. 
The two deep fields provide the only currently available sight- 
lines for HST data. We therefore need to use a model to obtain 
the color selection criteria that isolate high redshift star form- 
ing galaxies. 

The colors of galaxies as a function of age, metallicity, stel- 
lar initial mass function and star formation history can be ef- 
fectively modeled using the evolutionary population synthesis 
technique (see Bruzual & Chariot 2003, hereafter BC03, and 
references therein). Assumptions about the time evolution of 
these parameters allow one to compute the age dependent dis- 
tribution of stars in the Hertzsprung-Russell diagram, from 
which the integrated spectrum for the entire galaxy can be 
obtained. The models exploit the property that stellar popula- 
tions with an arbitrary star formation history can be thought 
of as the sum of a series of instantaneous starbursts (delta 
functions), referred to as simple stellar populations (SSPs). 
The spectral energy distribution at time t of a stellar popula- 
tion characterized by a star formation rate ip(t) and a metal- 
enrichment law ((t) can be written (following, e.g., Tinsley 
1980) 

F x (t)= [ i>(t-t')S x [t',at-t')] dt', (1) 
Jo 

where S\ [t'X(t-t')\ is the power radiated per unit wave- 
length per unit initial mass by an SSP of age t ' and metallicity 
C(f-f')- The above expression assumes that the initial mass 
function (IMF) does not change with time. 

The largest source of uncertainty in the spectral synthe- 
sis technique arises from the relatively poorly understood ad- 
vanced phases of stellar evolution, such as the supergiant and 
asymptotic giant branch (AGB) phases (BC03). 

In this paper, we use predictions from the BC03 model 
for the purpose of defining the color selection criteria for U 
dropout galaxies for the combination of ACS and WFPC2 fil- 
ters specific to our data. We chose to use the BC03 model, be- 
cause it is a recent model that provides high resolution spectra 
for our wavelength range of interest. The code implement- 
ing the model prescriptions, named GALAXYEV, is publicly 
available from the authors, thus providing a means of making 
our simulations reproducible by other researchers. 

We summarize here the main features of the BC03 model 
that are relevant for our purposes. The model predicts the 
spectral evolution of stellar populations of different metallic- 
ities and ages between 1 x 10 5 and 2 x 10 10 years at a reso- 
lution of 3 A FWHM over the whole wavelength range from 
3200-9500 A. These predictions are based on a new library of 
observed stellar spectra compiled by Le Borgne (2003) called 
STELIB. Predictions over a wider wavelength range but at 
lower resolution are also available. The Padova 1994 stel- 
lar evolution prescription and the BaSel 3.1 spectral calibra- 
tions are used. Two choices of IMF are available: the standard 
Salpeter (1955) IMF and the Chabrier (2003) IMF with lower 
and upper mass cutoffs of 0.1 M Q and 100 M Q respectively. 
Each model SSP is normalized to a total mass of 1 M Q in 
stars. The spectra are computed at 221 unequally spaced SSP 
ages from to 20 Gyr. Metallicity ranges from 0.005 to 2.5 
Z . 

The authors of BC03 have tested the predictions of their 
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FIG. 1 . — The final drizzled image in F300W obtained by combining 409 individual exposures. The white polygon encloses the area over which the wide field 
chips of all datasets overlap to give a total exposure time of 323.1 ksec, outside this area effective exposure is between 25% and 75% of the full depth. Portions 
of the image above and to the left of the black line lie outside the GOODS area and are not useful for identifying U dropout galaxies. 



model against observed color magnitude diagrams (CMDs) 
and integrated colors of star clusters. Overall, their model 
provides excellent fits to star clusters of different ages and 
metallicities in several photometric bands. They also com- 
pare the predictions of their model with observed galaxy spec- 
tra in the Early Data Release of the Sloan Digital Sky Sur- 
vey (SDSS). Comparison of the strengths of various absorp- 
tion features and Lick indices between the SDSS galaxies, the 
BC03 models, and other models shows a broad consistency 
between model and observations. 

The spectrum of a star forming galaxy, when observed over 
cosmological distances, is modified significantly by two ef- 
fects: (1) intergalactic attenuation along the line of sight from 
the galaxy to the observer and (2) dust attenuation in the emit- 
ting galaxy. In order to obtain a realistic model of the ob- 
served galaxy spectrum, it is necessary to account for both 
these effects. 

3.2. Intergalactic attenuation along the line of sight 

The model of Madau (1995) for the propagation of UV 
radiation through a clumpy universe is the standard work 
for determining intergalactic attenuation. This model for the 
H i opacity of the universe as a function of redshift includes 
scattering in resonant lines of the Lyman series (Lya, Ly/3, 
Ly7 and Ly£) and Lyman continuum absorption. Recently, 
Meiksin (2006) has proposed a model based on more current 
estimates of the properties of the intergalactic medium (IGM). 
The results of numerical simulations were used to estimate the 
contributions to resonant scattering from the higher-order Ly- 
man transitions. Differences of 0.5-1 mag from the previous 
estimate of Madau are found. In this work, we have chosen to 
use the standard work of Madau (1995) to enable us to make 



comparisons with previous work that have almost without ex- 
ception used this H i attenuation model. 

The amount of flux attenuation can be computed as a func- 
tion of wavelength at each redshift of interest. We show in 
Figure |2 a plot of the transmitted power fraction as function 
of wavelength and redshift. The characteristic staircase pro- 
file is due to continuum blanketing by the Lyman series. At 
z > 1.8, the flux in the f/-band is significantly lowered by 
the increased H i opacity. At higher redshifts, flux in progres- 
sively redder bands is affected. 

3.3. Attenuation by dust in the emitting galaxy 

We use the simple but realistic prescription of Chariot & 
Fall (2000) to account for the attenuation of light by dust. 
In this prescription, the attenuation of starlight by dust is ac- 
counted for by inserting a factor exp[-f>(f ')] in the integrand 
on the right-hand side of equation Q , where f\(t ') is the 'ef- 
fective absorption' curve describing the attenuation of pho- 
tons emitted in all directions by stars of age t' in a galaxy. 
This is given by the simple formula 

| f v (A/5500)^, for,'<10V, 
[ /if v (A/5500) , for t' > 10 7 yr, 

where fy is the total effective V-band optical depth seen by 
young stars. The characteristic age 10 7 yr corresponds to the 
typical lifetime of a giant molecular cloud. The adjustable 
parameter /i defines the fraction of the total dust absorption 
optical depth of the galaxy contributed by the diffuse inter- 
stellar medium (ji ps 1/3 on average, with substantial scatter). 

3.4. Our simulations 
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FIG. 2. — Mean transmission spectrum for a source at z=l. 8, 2. 8, 3.8,4. 8 (left to right, dotted lines) computed according to the model of Madau (1995). The 
characteristic staircase profile is due to continuum blanketing from the Lyman series. Also plotted are the response functions of our U,B,V filters (left to right, 
solid lines). 



The main aim in our simulations is to exploit the combined 
effect of the intrinsic Lyman edge in galaxy spectra and the 
opacity of intergalactic H i gas to separate relatively unob- 
scured star-forming galaxies from those that are older and/or 
dustier, in a color-color diagram of available broad-band col- 
ors. Such a separation in color-color space is most effectively 
achieved by plotting the U-B color against the B-V or B-I 
color. In practice, the choice of B- V or B-I does not have 
a significant impact on the efficiency or completeness of the 
separation. Following COO, we chose the B-V color as the 
redder color in our simulations. 

We wish to compute colors for galaxies populating a grid of 
ages, metallicities, dust obscuration and star formation histo- 
ries at various redshifts (see Table^for a summary of param- 
eter grid values used in our simulations). As a first step, we 
need to study how the choice of IMF (Salpeter or Chabrier) af- 
fects the distribution of galaxies in the color color space. The 
Chabrier IMF is identical to the Salpeter IMF for stars that are 
more massive than 1 M Q . For lower mass stars and substellar 
populations such as brown dwarfs the two IMFs are signifi- 
cantly different. However, the rest-frame UV emission from 
galaxies at high redshift that is probed by the colors used in 
our simulations, receives a nearly negligible contribution from 
low mass stars and brown dwarfs. Thus, for our simulations, 
choosing either IMF should produce nearly indistinguishable 
color-color plots. We verified this by comparing color color 
plots for both IMFs in the relevant region of color-color space. 
For consistency with previous work, we chose to use the stan- 
dard Salpeter IMF in our models. There is no prior on which 
points on the grid are valid, except the obvious logical re- 
quirement that the galaxy should be younger than the age of 



the universe at that redshift. The grid thus samples the plausi- 
ble range of galaxy properties, without weighting the results 
with our knowledge of the actual distribution of galaxy color 
that may be derived from either observations or theoretical 
models. 

We compute the spectral energy distribution (SED) of each 
galaxy by modeling it as a composite stellar population that 
formed either in a constant burst of star formation lasting 10 7 
years or in an exponentially decaying burst with timescales 
ranging from 10 6 years to 3 x 10 10 years. The SEDs are com- 
puted at a range of ages (defined as the time since the onset of 
star formation) and with a range of metallicities. Note that all 
stars in a given galaxy have exactly the same metallicity in the 
BC03 model. Due to this, the change in integrated colors due 
to metallicity gradients within a galaxy are not model ed. The 
SEDs are corrected for dust attenuation (Section f3.3t . appro- 
priately redshifted and then flu x cor rected for cosmic opacity 
along the line of sight (Section f3.2t . Finally the SEDs in the 
observer frame are folded through the filter response func- 
tions and the broad-band colors in the U-B and B-V filter 
are obtained. 

3.5. Colors of relatively unobscured star forming galaxies at 
2.0<z<3.5 ' 

We show in Figure [3] the color-color diagram that we ob- 
tained for our simulation of galaxies using the parameter grid 
described above. Colored symbols (cyan through magenta) 
identify galaxies within the desired redshift range 2 < z < 3.5 
that are both starforming (age < 0.1 Gyr) and relatively un- 
obscured by dust (T\(t') < 2.0). Brown symbols identify old 
and/or attenuated galaxies in the same redshift range; such 
galaxies are typically much fainter in restframe UV (observed 
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TABLE 1 




Parameter grid used for computing synthetic spectra. 




Ages 


0.001, 0.0025, 0.005, 0.0076, 0.01, 0.025, 0.05, 0.1, 0.5, 1.0, 5.0, 10.0 Gyr 


Exponential star formation timescales 


0.01,0.1, 1.0, 5.0, 10.0, 30 Gyr 


Constant star formation durations 


0.01 Gyr 


Metallicities 


0.2, 0.4, 1.0 Z Q 


rx{t') 


0.01,0.51,1.01,1.51,2.01,2.51 


M 


0.02, 0.22, 0.42, 0.62, 0.82 


Redshifts 


0.2-7.0 Stepsize: 0.2 




■1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

B-V 

FIG. 3. — U -B versus B — V for our simulated galaxies. A total of 264 600 synthetic spectra of galaxies representing a wide range of ages, star formation 
histories, metallicities, dust properties and redshifts were folded through our combination of WFPC2 and ACS filters. Galaxies shown in colors other than black 
represent young star-forming galaxies with relatively low level of dust obscuration in the redshift range 2.0 < z < 3.5. The colors cyan, blue, green, orange, 
red and magenta show galaxies with progressively higher redshift in bins of 0.25 per color. Galaxies in the same redshift range that are old and/or severely dust 
attenuated are shown as gray dots. Black dots represent galaxies outside the redshift range. Our color selection criterion is shown by the polygon. Galaxies 
within the polygon are the U-band dropout candidate galaxies with 2.0 < z < 3.5. The selection criteria are U — B > 1.0, U—B > B- V + 1.3, and B — V < 1.2. 
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B) and are unlikely to constitute a large fraction of our sample. 
Black dots identify galaxies outside the target redshift range. 

We find that relatively unobscured star forming galaxies at 
2 < z < 3.5 are effectively selected in color-color space by the 
following color criteria: (1) U-B > 1.0 (2) U-B > B-V + 
1.3 and (3) B-V < 1.2. The minimum "redness" threshold 
in U - B color (horizontal line) is motivated by the desire to 
exclude objects with z < 2. For a given metallicity and dust 
attenuation, galaxies with 2 < z < 3.5 become redder in U — 
B at progressively higher redshifts while remaining relatively 
unaffected in their blue B-V color. As the redshift approaches 
z < 3.5, the Lyman a forest and eventually the Lyman limit 
begins to progressively intrude into the B band causing the 
B-V color to become redder. This begins to move galaxies 
rightward and out of the selection box. Increasing metallicity 
and dust obscuration both tend to make objects redder; the 
extent of the effect on each color is strongly dependent on the 
age of population and the redshift. 

By design, these selection criteria are aimed at identifying 
galaxies with strong UV flux. Therefore our sample will ex- 
clude galaxies with a low star formation rate or with strong 
reddening. The galaxies to the right of our B-V cutoff in 
Figure [3] are relatively dusty {f\{t') > 1.0) and are forming 
stars at a relatively low rate. Such moderately reddened galax- 
ies with modest star-formation rates have been identified and 
studied spectroscopically at brighter magnitudes at somewhat 
lower redshift (Daddi et al. 2004; Reddy et al. 2005), but their 
luminosity function has not been probed directly to the faint 
levels we are investigating here. Similarly, recent ground- 
based studies (e.g., Le Fevre et al. 2005) find a non-negligible 
number of spectroscopically confirmed high redshift galaxies 
below the selection box (bluer U-B color). Again, these ob- 
jects are significantly brighter (/ < 24) than the objects that 
we aim to study in this work. Significantly expanding the 
selection box in either direction would include many more 
contaminating objects outside our target redshift range; there- 
fore we prefer to maintain a conservative selection criterion 
that keeps contamination low, while including the majority of 
strongly star-forming, relatively unobscured galaxies. This 
choice is consistent with that of other dropout studies, al- 
though our color limits are slightly different from COO due 
to the updated stellar population models and to the fact that 
the ACS B (F435W) filter we use is shifted blueward of the 
WFPC2 B (F450W) filter. 

3.6. High redshift galaxies in the HUDF parallels: The 
sample 

Having identified the color selection criteria above, we are 
now in a position to identify U dropout galaxies in our data. 
To obtain photometry for our images, we used Sextractor ver- 
sion 2.3.2 (Berlin & Arnouts 1996) in dual image mode. We 
used the GOODS/ACS i band image for source detection and 
the WFPC2 £/-band and GOODS/ACS B and V band images 
for photometry. The region of the GOODS data that fell out- 
side the footprint of the U band image was set to zero. 5.67 sq. 
arcmin (87%) of the 6.51 sq. arcmin footprint of our [/-band 
image overlaps with the GOODS coverage area (see Figure^. 

The main advantage of dual mode photometry is that there 
is no need to make a subjective decision to match correspond- 
ing sources detected in the various bands. Setting up such a 
correspondence is particularly difficult in situations where a 
single object in one filter gets detected as multiple objects in 
another filter, if there are multiple peaks in the light distri- 
bution. Such a situation is common in the U band which se- 



lectively probes the regions of ongoing star formation within a 
galaxy. Detection was carried out with a DETECT_MINAREA 
of 12 pixels and a DETECT_THRESH of 1.5. The detections 
were filtered with a Gaussian kernel with FWHM 2.5 pix- 
els and cleaned with CLE AN_P ARAM= 1.0. Object deblending 
was carried out with DEBLEND_MINCONT=0.03 

The sample is first truncated at magnitude limits sufficiently 
bright that Lyman breaks of the expected amplitude could be 
measured reliably. For our data, this corresponds to B < 27.3. 
We find that all objects that are brighter than this limit in 
the B band are also detected in the i band image. We show 
in Figure |4] the color-color plot highlighting the candidate 
[/-band dropouts. Objects with measured magnitudes in all 
three bands are indicated by crosses. Objects with measured 
U band flux (as indicated by the Sextractor MAG_AUTO mag- 
nitude) fainter than the la limit computed from the inverse 
variance map at that location are considered undetected, and 
are marked by a triangle at the U-B color corresponding to 
the la magnitude limit. The la limit depends on the size of 
the object. Additionally, as the exposure time is not uniform 
across the image, it also varies with position. The size of each 
symbol is proportional to the ; band magnitude of the object. 

The polygonal selection regio n wh ich identifies U band 
dropouts, as defined in Section 13.51 contains 136 objects, 
each of which has been visually inspected. Five objects, 
all relatively bright (i < 22), appear to be stars; 6 more ob- 
jects are otherwise spurious, either due to instrumental arti- 
facts (diffraction spikes, image edges) or resolved segments 
of nearby bright galaxies. These 1 1 objects are not plotted in 
Figure |4] and have been excluded from the subsequent anal- 
ysis. The remaining clean sample of U dropouts consists of 
125 objects; their coordinates and magnitudes are listed in Ta- 
ble|2] 

4. DROPOUT MORPHOLOGY 

At redshifts z ~ 3, the drizzled ACS images with a resam- 
pled pixel scale of 0.03 arcsec/pixel probe regions as small 
as 0.23 kpc/pixel for our chosen cosmology. The fine resolu- 
tion provides sufficient resolution elements to enable at least a 
crude study of their morphology. The drizzled GOODS/ACS 
images have the advantage of a pixel scale that is a factor 
of two finer than the drizzled HDF images. However, the 
depth reached by the GOODS survey is about a magnitude 
shallower than the corresponding HDF data. Morphological 
details can thus be studied with better resolution but poorer 
signal-to-noise, relative to the HDF images. 

We have visually inspected the entire sample of U -dropouts 
in the original GOODS/ACS z-band images and classified 
their morphology using an empirical classification scheme. 
Objects that appear to be single and compact are classified as 
clumps (1), double or multiple clumps (chains) are classified 
as (2), clumps with tails, i.e. tadpoles as (3). Objects which 
are not compact but extended and fuzzy are classified as (4), 
and fuzzy extended objects with a more organized disk-like 
structure are classified as (5). Figure [5] illustrates typical ob- 
jects of each class. We list in Table[3]the relative fraction of 
each class. For our sample, the fraction of single clumps is 
similar to the fraction of fuzzy objects (~ 35%), however if 
double clumps and tadpoles are included in the clump class 
(1), this fraction increases to ^50%. The morphological clas- 
sification for each object and short remarks on its appearance 
are listed in Table |2] We also show in Figure |6] BVi color 
composites of 16 typical, bright dropouts in our sample. 

Since the original deep fields it has been known that com- 
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FIG. 4. — Color-plot indicating the candidate (/-band dropouts (enclosed by the polygon) in our data. Crosses indicate detections; triangles indicate non- 
detections in the J/-band and are placed at the la lower limit of the U-B color. The size of the symbols scales with magnitude in ACS/F775W. 



pact/clump objects dominate at faint magnitudes (Abraham et 
al. 1996, van den Bergh 2002). Recently, this result has been 
confirmed by the morphological studies of Ultra Deep Field 
where most of the galaxies are found to be clumps (Elmegreen 
et al. 2005a, 2005b). An important question is whether such 
clumpy structures are signatures of merging/accretion or a 
merely band-shifting effect. The latter is probably not the 
case since these clumps have no counterparts in the local Uni- 
verse (Elmegreen et al. 2005a, 2005b, Windhorst et al. 2002). 
Moreover, as shown in Papovich et al. (2005) for galaxies at 
1.9 < z < 3, there is very little morphological transformation 
between wavelengths using WFPC2 F606W image and NIC- 
MOS F160W which spans a long wavelength baseline (rest- 
frame far-UV to the rest-frame B band). 

We have also checked whether the morphological counter- 
parts of U -dropouts are seen in the lower redshift (z < 2) sam- 
ple (de Mello et al. 2006) for which we have photometric red- 
shifts from GOODS (Dahlen et al. 2005). We found that U- 
dropout objects show morphological similarities to ones clas- 



sified in de Mello et al. (2006) as starbursts 5 which are clas- 
sified as 'compact, peculiar or low surface brightness (LSB)'. 
The objects classified as 'compacts' are similar to class (1). 
Objects classified as (2) and (3) would fit into the de Mello et 
al. 'peculiar' class, and objects classified as class (4) (fuzzy) 
correspond to the LSB class. However, as expected, the U- 
dropout sample lacks the classical Hubble types, such as el- 
liptical/spheroid and disks which are seen at lower redshifts 
(z < 0.8). The fraction of LSBs and compact+peculiar ob- 
jects in the sample of de Mello et al. (2006) which are at 
0.8 < z < 1.2 is 27% and 52%, respectively. Interestingly, the 
corresponding fractions of [/-dropouts are very similar, 38% 
and 51% (extended fuzzy and clump+doubles+tadpoles) in- 
dicating that about 50% of the star formation takes place in 
clumps over an extended redshift range from 0.8 < z < 3.5. 



5. DROPOUT STATISTICS 

5 Spectral type obtained from the template fitting in the photometric red- 
shift technique (Dahlen et al. 2005). 
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FIG. 5. — Representative dropout in each of the 5 classes into which we classify objects. The classes are (1) clumps, (2) double or multiple clumps, (3) clumps 
with tails, i.e. tadpoles, (4) objects which are not compact but extended and fuzzy and (5) fuzzy extended objects but with a more organized disk-like structure. 



It is interesting to compare the statistical properties of the 
dropouts selected here with those in the two HDFs. The selec- 
tion procedures are very similar, although slight differences 
can arise from the different filter set available and from the 
different depth of the observations — those used here are sig- 
nificantly deeper in U, but shallower in the other bands. With 
those caveats in mind, we compute some properties of interest 
for our sample. 

5.1. Number density 

We have identified 125 U band dropouts in the 5.67 square 
arcmin of sky over which we have multiwavelength data avail- 
able. Of these, 80 fall within the full depth coverage of our 
F300W data which covers 4.02 square arcmin. The surface 
density of dropouts at the full depth of our study is thus 19.9 
sources per square arcmin, similar — within the statistical 
uncertainties — to the 17.5 sources per square arcmin found 
in the HDF-S data (COO). Relative to the HDF-N value of 
16.0 sources per square arcmin our measured surface den- 



sity is higher with a marginal significance. Both the HDF 
measurements were obtained using the Madau et al. (1996) 
color selection criteria and only included sources brighter than 
B450 = 26.79. Our sample uses a different filter combination 
with its own color selection criterion and all sources brighter 
than B = 27.3. 

Using only the dropouts in the region with full depth and as- 
suming that these dropouts uniformly probe the redshift range 
2 < z < 3.5, we obtain a comoving galaxy density of 4.0 x 
10~ 3 Mpc~ 3 for the concordance cosmology we adopt, cor- 
responding to 5.7 x 10~ 3 Mpc" 3 with the cosmology adopted 
by Madau etal. (1996; qo = 0.5,/iioo = 0.5, no cosmological 
constant). For comparison, they report a density of 4.2 x 10~ 3 
Mpc" 3 from the HDF-N. Our somewhat higher galaxy den- 
sity may be due in part to the increased depth of our sample, 
although cosmic variance can also play an important role on 
such small areas. 

5.2. Star formation rate density evolution 
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FIG. 6. — BVi color composite of 16 bright U band dropouts in our sample. Each cutout is 5x5 arcsec in size. The dropout galaxy is located at the center of 
each panel. 



A key measure is the evolution of the cosmic star forma- 
tion rate (SFR) density as a function of lookback time (or 
redshift) (e.g., Madau etal. 1996). The star formation rate 
may be measured using carefully calibrated proxies such as 
the luminosity in X-ray (Ranalli et al. (2003), radio (Condon 
et al. 1992), far-infrared (Kennicutt 1998), narrow-band line 
emission (Moustakas et al. 2006) and the rest frame UV con- 
tinuum (Madau et al. 1998). When such indicators are used 
to compute the SFR for the same sample, they are consistent 
to within a factor of two (Daddi et al. 2004). 

Even within various samples identified using the dropout 
technique, significant discrepancies may exist. These may be 
caused by cosmic variance, differences in stellar population 
synthesis models or in the treatment of dust attenuation. Dif- 
ferent assumptions on the shape of luminosity function and 
the differing photometric uncertainties in observations of dif- 
ferent depths can also introduce systematics into the mea- 
surement. Of these, the uncertainties in the dust obscuration 
correction are of the greatest concern in estimating the SFR 
from the rest frame UV continuum. Recently, Heckman et 



al. (2005) have obtained GALEX observations of compact 
UV luminous galaxies in the local universe. They find that 
these galaxies, which are chosen to have rest-frame UV lu- 
minosities similar to those of z ~ 3 star forming galaxies also 
match them in stellar mass, velocity dispersions and gas phase 
metallicity. Like their distant counterparts, these galaxies are 
forming stars at a rate high enough to build their present day 
stellar mass in 1-2 Gyr. For these local galaxies, dust attenu- 
ation is relatively modest - 0.5-2 magnitudes. If the high red- 
shift galaxies exhibit similar dust properties, the effect of un- 
certainty in dust attenuation on the computed star-formation 
rate may be less than previously thought ( Vijh, Witt & Gordon 
2003). 

In spite of these uncertainties, general trends in the evolu- 
tion of the SFR density as computed from the UV continuum 
luminosity have become apparent over the last decade. We 
now know with a high level of confidence that the SFR den- 
sity steadily increases from the local Universe out to a redshift 
of z ~ 3 and then either shows a fall off or remains steady 
at higher redshifts (e.g., Lilly et al. 1996; Madau et al. 1996; 
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FIG. 7. — SFR density evolution as a function of redshift (compare to Figure 11 of COO). The SFR is estimated from the UV luminosity density. The z > 2 
points are from Lyman-break objects in the HDF-N (open triangles), in the HDF-S (filled triangles) from COO, and in the Steidel etal. (1999) ground-based 
survey (x symbols). The black ellipse indicates the SFR density measured in this work. For clarity, error bars are not plotted for our datapoint. Errors along 
the redshift axis are identical to the HDF points. The statistical error on the SFR density axis is about 0.09. Distances and volumes are computed using the 
cosmological parameters h, f! ra , Ha, f2tot = 0.65,0.3,0.7, 1.0.. At lower-redshifts, only results that are directly comparable to the high redshift data are plotted. 
The open squares are from HDF photometric redshifts by Connolly et al (1997), the open circles are from Lilly etal. (1996), the solid square is from the Ha 
survey of Gallego et al. (1995), and the solid circle from Sullivan et al. (2000). 



Dickinson 1998; Steidel et al. 2003; Bouwens et al. 2004). We 
adopt here the formulation of COO, which in turn follows the 
prescriptions of Steidel et al. (1999). For ease of comparison 
with Figure 1 1 of COO, we use their cosmological parameters 
/z,f2 m ,f2A,^tot = 0.65,0.3,0.7, 1 and define a fiducial magni- 
tude /?606+775 for each candidate, which is the average of the V 
and i magnitudes. We use this magnitude to compute n?AB,tot 
which is the total integrated flux of the Lyman-break candi- 
dates. To avoid variable levels of incompleteness, we only 
use dropouts in the full depth region of our image for comput- 
ing the integrated flux. Then, the star-formation rate density 
at a mean redshift z between two redshifts zi , Zi is given by 



log( 



SFR 



M Q yr _1 Mpc~ 



f ) = 21og 
-log 



D L (z) 



-0.4mA 



cm 

AV( Zl ,z 2 ) 



log A 



Mpc 3 

+log<SL-log(l+z)-34.516, 



(3) 



where the redshifts z,zl,z2 are 2.75,2.00,3.50respectively, Dl 
is the luminosity distance, AV is the total volume between 
the two redshifts, and A is the area of the field of view in 
square arcseconds. The contribution of galaxies fainter than 
our detection threshold is included via the correction factor 
6L, which is obtained by integrating the assumed luminosity 
function (Steidel et al. 1999) from our survey limit down to 



0.1L* (Formula 9, COO). The above formula includes a uni- 
form extinction correction factor of 4.7 (Steidel et al. 1999) to 
the measured UV luminosity densities. This correction factor 
is appropriate for a typical extinction E(B-V) = 0.15 and the 
Calzetti (1997) reddening law. 

We measure an integrated star-formation rate density of 
0. 18 M0yr _1 Mpc -3 . This value is only marginally higher than 
that measured for the HDF-S, which in turn is higher than that 
measured for the HDF-N. Both these values were measured 
by the authors of COO, using identical source detection and 
photometry procedures with data of similar depth obtained 
using the same set of filters and with the same color selection 
criteria. We show in FigureQthe evolution of star formation 
rate density with redshift, comparing our measurement with 
other measurements from the literature obtained using a sim- 
ilar technique. 

Measurements of the statistical properties of any clustered 
population are plagued by uncertainty due to cosmic vari- 
ance, the field to field variation due to large scale structure 
(Somerville et al. 2004). At high redshifts, especially with 
the small few arcmin sized field of view of HST, the volume 
sampled between 2.0 < z < 3.5 are ~ 10 4 Mpc 3 , small enough 
that cosmic variance is a significant source of uncertainty. The 
variance of number counts N over an area A on the sky, in 
the presence of a non-zero angular auto-correlation function 
is given by (Peebles, 1980): 
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a 2 (N) = NA + n 2 



dA x / dA 2 w(6) 



(4) 



where n is the average density of objects and thus nA is the ex- 
pectation value of N, and w(9) is the two-point angular auto- 
correlation function between points A\ and A2. For our full 
depth area of 4.02 sq. arcmin the second term in equation [5] 
can be approximated as \.%n 2 A 2 A w where A w is the two point 
angular correlation of the sample at a scale of 1 arcmin. This 
assumes that the angular correlation function obeys a power 
law with exponent [3 = — 1.2 as measured by Giavalisco & 
Dickinson (2001) for a combined sample of U dropouts from 
the two HDFs. We estimate the value of A w at 1 arcmin from 
the la upper limit on their measured angular correlation func- 
tion for their combined HDF sample as A w = 0.02. The total 
expected variance is thus: 



a 2 (N) = 80+ 1.8 x 0.02 x80 2 = 310.4 



(5) 



or approximately 3.9 times larger than the Poisson contribu- 
tion alone. The expected r.m.s. fluctuation in number counts 
is about 17.6 galaxies over the full depth area, compared with 
8.9 for the Poisson case. This corresponds to an uncertainty 
of 0.09 in the logarithm of the star formation rate density, and 
may be treated as the nominal statistical uncertainty on our 
measurement. It should be noted that this error bar only ac- 
counts for statistical uncertainties; systematic offsets e.g. due 



to incorrect dust corrections may be significantly larger. 

An effective way to reduce the uncertainty caused by clus- 
tered populations is to obtain samples of high redshift galax- 
ies along independent lines of sight. Our measurement of the 
SFR density in the GOODS-S area provides such an inde- 
pendent measurement, and is in fact only the third such mea- 
surement obtained exclusively from space-based data at z ~ 3. 
Within the error bars, our measurement of SFR density is con- 
sistent with the value measured for HDF-S, indicating that the 
range of variation in SFR density is probably well sampled by 
the two HDFs. 

In the future, the second WFPC2/F300W parallels, if com- 
bined with additional new data in the B,V filter can provide a 
fourth space-based measurement of the SFR density at z ~ 3, 
further reducing the impact of cosmic variance. 
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9A 9A 


93 QQ 

zj.yy 


9 

z 


double clump 


o 171 17/1 
JJ.l / 1 J /4 


97 70QAAA 

-2.1 . /uy44o 


98 89 
Z5.5Z 


9A AA 
Z0.44 


9A OQ 

zo.uy 


9S 03 

Zj.yj 


1 
J 


faint clump 


Si i sasso 

j j. 1 jOjjU 


97 70S 81 7 
-Z / . /U5o 1 / 


98 A8 
Zo.^fo 


9A OA 

zo.yo 


9A 77 
ZO. / / 


9 A SO 
ZO.OU 


1 

1 


faint clump 


Si i AA9AA 

Jj. lOOZOO 


97 7n8SAS 
-Z / . /UojOj 


98 08 
zo.Uo 


9A 93 
ZO.Zj 


9A OA 
ZO.UO 


9A 08 
ZO.Uo 


1 

1 


clump 


Si i SS7A9 
JJ. 1 JJ /oz 


97 70R30S 
-Z / . /UojUj 


90 S3 

zy. j j 


9A A7 
ZO.4- / 


9S AS 
ZJ.4-0 


9S 1 9 
Zj. 1Z 


s 

j 


clumps in disk 


O 1 AAA1 1 
Jj.l 0O0 1 1 


97 708980 
-2.1 . /U5Z5U 


98 S7 
Zo. J / 


9A S3 
ZD. J j 


9A IS 
Z0. Jo 


9S SQ 

zj.sy 


4 


faint fuzzy 


SI 1 A0SS7 
JJ. lOUJJ / 


97 707AAQ 

-Li . ivj i ooy 


97 3S 

Z / .JJ 


9S AS 
Zj.Oo 


9S 33 
ZJ. jj 


9A 77 
Z^f . / / 


9 
Z 


multiple clump (maybe on faint disk) 


Si i 10A9A 

JJ. IjUDZH 


97 707Q1 9 

-z/. /u/yiz 


>, 98 SA 
S> Z5.5D 


97 1 9 
Z / . 1Z 


97 00 
Z / .UU 


9A 3Q 

zo. jy 


s 

J 


faint clumpy disk/tail? 


si i soaoq 
j j. i juouy 


97 707880 
-Z /. /U/ooU 


97 on 
z / .yu 


9^ A7 
Zj.O / 


9S /IS 
ZJ. 45 


9^ An 

Zj.4 / 


A 

4 


faint fuzzy 


Si i S3A8A 

J J. 1 J J054 


97 7071 SA 
-Z / . IV) 1 1 50 


97 3A 
Z / . j4 


9S 9S 
Z J.Zj 


9A 83 
Z4-.0 J 


9A A3 
Z^+.Oj 


A 
4 


f uz zy/clump 


S3 1 Sn37A 
Jj. IjUj /4 


97 70A89A 
-Z / . /UOoZD 


97 SO 

z / . jy 


9S OS 

Zj.Uj 


9A 7A 
Z4-. /'f 


9A AA 
Z^f .^fO 


s 

j 


faint clumpy disk 


S3 i n OH3 
JJ. IjIUUj 


97 70AAQ8 
-Z / . /U04yo 


9A OA 

zo.y4 


9S 97 
Z J.Z / 


9S 93 
Z J.Z j 


9A AS 
Z^f .OJ 


1 
1 


clump 


S3 1A30QA 
j j. i4juyo 


97 70SQSA 

-z / . /ujyjo 


90 AQ 

zy.oy 


9 A 80 
ZO.OU 


9A SA 
ZO. J^f 


9A SO 
ZO. JU 


s 

J 


faint clumpy disk 


S3 17SS7Q 

j j. i / jj / y 


97 70S789 
-Z / . / Uj / oZ 


97 3S 

Z / .JJ 


9S 90 
Z J.ZU 


9A 89 
Z4-.OZ 


9A A9 
Z'+.OZ 


9 

z 


double clump - merging 


S3 i inns a 

J j. 1 jUUOD 


97 70S7S8 
-Z / . / UJ /Jo 


98 80 

^> Zo.oy 


9A AQ 

zo.oy 


9 A Al 
ZO.4-1 


9A 3 1 
ZO.J 1 


A 
4 


faint fuzzy 


1 711 in 

jj.l / j 1 jU 


97 70SA3 1 
-Z / . /UjDj 1 


98 01 

zo.y i 


9A 77 
ZO. / / 


9A 98 
ZO.Zo 


9 A 07 
ZO.U / 


A 


faint fuzzy 


SI 1 89801 
Jj . 1 oZoU 1 


97 70S979 
-Z / . /UjZ / Z 


97 no 
z / .uy 


9A S7 
Z4-. J / 


9A S3 
Z4-. J J 


9A SS 
Z^f . J J 


1 

1 


clump edge on 


SI 1 1SQ98 
jj. i jjyZo 


97 70S9QQ 

-z / . /ujzyy 


90 00 

zy . uu 


9A Q7 

zo.y / 


9A QS 

zo.yo 


9A SA 
Z0.50 


A 
4 


faint fuzzy 


SI 1 8 1 QQ1 

JJ . 1 O 1 77 1 


97 70S09Q 
-z / . /ujuzy 


28 49 


26 90 


26 72 


26 76 


4 




S3 150622 

JJ.l JUUiL 


-27 704625 


27.38 


25.40 


25.35 


25.32 


I 


clumn nTTPfntfT" 


SI 17/1787 
JJ.l /4/0 / 


97 70AAQQ 

-Z /. /U44yy 


10 1 1 
> jU. 1 1 


9A AS 
ZO.Oj 


9 A 09 
ZO.UZ 


9 A OQ 

zo.uy 


A 
4 


faint fuzzy close to spiral 


SI 1AS378 
JJ. ItOJ / o 


97 70A9Q8 
-Z / . /U4Zyo 


90 81 

zy.o i 


9A OA 
ZO.UO 


9S 09 
ZJ.UZ 


9A 79 
Z^f . / Z 


1 
1 


clump (spheroid) 


SI 1 ASAQ8 
JJ. IDJD70 


97 703QS0 

-z / . /ujyju 


\ 9Q Q0 


9A 83 
ZO.OJ 


9A 1A 
ZO. l*f 


9S QA 
ZJ.y4- 


A 
4 


very faint fuzzy 


SI 111 1 AQ 

jj. i j i loy 


97 703QA8 
-Z / . /UjyOo 


97 SA 
Z / .oO 


9S as 

ZJ.OJ 


9S 80 
Zj. oU 


9S 7A 
Zj. Ih 


1 
1 


clump/edge on 


SI 1 10S01 
JJ. IjUjUI 


97 703 AGO 
-Z / . /ujoyu 


9S 37 
Zo.j / 


9A 93 
ZO.Zj 


9 A AO 
ZD.4U 


9S 77 
Zj. / / 


s 

J 


fuzzy disk 


SI 17SQQ1 
Jj.l I jyy 1 


97 703 1 3S 
-Z / . / Uj 1 j j 


>, io A3 
^> jU.4j 


9 A Ql 

zo.y 1 


9A AA 
ZO.OO 


9A A7 
ZO.O / 


1 
1 


faint clump 


SI 140089 
jj. i4yu5Z 


97 709QS9 

-z / . /uzyjz 


98 7S 
Zo. / J 


9 A 09 

zo.uz 


9 A 90 
ZO.ZU 


9A 9A 
ZO.Z4- 


1 

1 


clump/edge on 


SI 1 A1 871 
jj. 141o / 1 


97 709771 
-Z / . /UZ / / 1 


98 AA 
Zo.^t 


97 1 1 
Z / . 1 1 


97 31 
Z/.jl 


97 01 
Z / . U 1 


A 
4 


fuzzy 


SI 1 377SA 
JJ. 1 J / / JO 


97 701 9Q9 

-z / . /uizyz 


>. 90 AS 

j> zy.oj 


9A 9A 
ZO.ZO 


9S as 

ZJ.OJ 


9S 71 
Zj. / 1 


1 
1 


faint clump offcenter on fuzzy tail (companions) 


SI 1 90 SAO 

jj. izyj4u 


97 701 970 

-z / . / uiz / y 


98 OA 

j> zo.yo 


97 93 
Z / .Zj 


97 00 

z / .uy 


9A 7A 
ZO. /*+ 


1 
1 


clump next to peculiar shape disk (multiple clump?) 


SI 1/11/117 
jj. 1414J / 


97 701 1 SO 
-LI. /Ul lJ7 


98 00 

zo.yy 


97 1 1 

Z / . 1 1 


9A QQ 

zo.yy 


9A S3 
ZO. J J 


A 
O 


very faint 


SI 1 3SSA1 
JJ. 1 jooO 1 


97 70 1 081 
-Z / . /UlUoj 


97 9A 
Z / .ZO 


9S A7 
Zj.O / 


9S Al 
Zj.41 


9S m 

Zj.UI 


1 


clump 


S3 178AS8 

JJ.l/ OtJO 


97 701 1 73 
-z / . / U 1 1 / J 


28 77 


27 23 


97 08 

Z / .UO 


26 33 


5 


faint disk edge-on 


SI 1 3S01 9 
J J. 1 JoUlZ 


97 70 1 1 91 
-Z / . /Ul IZj 


^ 90 01 

> zy.yi 


9A ^ 
ZO. JJ 


9 A 07 
ZO.U / 


9^ Q1 

zj.y 1 


1 
1 


clump/close to other faint clumps 


S3 1£1 S1 S 

JJ. lUlJlJ 


97 A8AQ7Q 
-z / .oo^+y / y 


98 no 

Zo.UU 


9A OA 

ZO.UO 


25 77 


25 59 


| 


clump/tail close to fuzzy 


SI 1 Al iss 

JJ. lOl JJ J 


97 A8S 1 1 A 
-Z / .OoJ 1 lO 


97 33 
Z / .jj 


9S AS 
ZJ.OJ 


9S AO 

zj.oy 


9S SA 
Zj. 50 


A 
4 


fuzzy close to clump/tail 


SI 1ASA30 
JJ. IIjDjU 


97 A8S9A1 
-Z / .OoJZOl 


97 oo 
S> Z / .oo 


9S Al 
ZJ.^tl 


9A AS 
Z4-.OJ 


9A AA 


1 
1 


c lump( spheroid) ? 


SI 1 ASA 1 Q 

jj. ioj4iy 


97 A8SA97 
-Z / .OoJOZ / 


98 87 
Zo.o / 


9A SQ 

zo. jy 


9A 1 

zo. iy 


9 A 07 
ZO.U / 


A 
4 


fuzzy 


SI 1A9Q79 
jj. 14Zy 1 Z 


97 A8SAQ8 
-Z / .OoJOyo 


9A 03 

zo.yj 


9S A7 
ZJ.4- / 


9S AQ 

zj.4-y 


9S 97 
Zj.Z / 


s 

J 


faint disk edge-on 


SI 1A0890 
JJ. 14UOZU 


97 A8S880 
-Z / .OoJooU 


^ 90 no 
j> zo.uy 


9 A AO 
ZO.OU 


9A AS 
ZO.OJ 


9A A7 
ZO.^f / 


s 

J 


fuzzy/ disky 


SI 1A3A0A 
JJ. 14 J4U4 


97 A8A93 1 
-Z / .OoOZj 1 


\ 98 79 
J> Zo. / z 


9A 83 
ZO.oJ 


9A AS 
ZO.Oo 


9A S3 
ZO. J J 


9 
Z 


double clump 


SI 1A7Q07 

jj. 14 /yu/ 


97 A8770A 
-Z / .Oo / /U4 


97 OS 

z 1 .yj 


9A 97 
ZO.Z / 


9A 1 S 
Z0. 1 J 


9A 9A 
ZO.ZO 


9 

z 


big chain 


SI 1AA981 
JJ. 140Z5 1 


97 A883A3 
-Z / .OoojOj 


\ 90 Al 
S> Zo.^f/ 


9A A8 
ZO.Oo 


9A A7 
ZO.4- / 


9A 33 
ZO. j J 


A 
4 


fuzzy close to fuzzy objects 


SI 1A78SS 
JJ. 14- / oJJ 


97 A88SQ1 

-z / .Ooo jy i 


98 37 
S> Zo.j / 


9A AA 
ZO.^t^f 


9A 39 
ZO. JZ 


9A 91 
ZO.Z 1 


9 
Z 


chain 


S3 1 A1 3AA 
JJ. 101 J40 


97 A8883A 
-Z / .Oo5o JO 


98 73 
Zo. / J 


9 A 30 
ZO. jU 


9S 83 
ZJ.OJ 


9S SS 
Zj.jj 


A 
4 


fuzzy /edge on 


SI 1 SAA9A 
J J. 1 J40Z4 


97 A8Q/1 1 8 
-Z / .Ooy41o 


\ 3n 11 

> JXJ.J 1 


97 1 A 

Z / . 10 


97 99 
Z / .ZZ 


97 IO 
Z / . JU 


1 
1 


faint clump 


SI 1A30AQ 
JJ. 14jU4y 


97 A8Q79A 

-z / .ooy /Z4 


>, 98 A9 
J> Zo.OZ 


9A 38 
ZO. Jo 


9A 9S 
ZO.ZJ 


9A 9 1 
ZO.Z 1 


A 
4 


fuzzy 


SI 1 AQ9Q9 

jj. loyzyz 


97 A8QQ70 

-z / .ooyy /u 


98 38 
Zo. jo 


9A OA 
ZO.U4 


9S S3 

Zj.jj 


9S AO 
Zj. 'tU 


1 
1 


clump 


si i/i^/i no 

jj. 14J4Uy 


97 A00901 

-z / .oyuzuj 


97 Q/1 

z / .y4 


9S 8Q 

zj.oy 


9S /IS 
ZJ. 45 


9S /IO 
Zj.4U 


9 
Z 


clump 


SI 1AS39S 
JJ. IH-jjZj 


97 AQ09S9 

-z / .oyuzjz 


29 20 


26 06 


25 47 


25 37 


2 


clump 


si 1 ^7^^n 

JJ. 1 J / jjU 


97 Aon/ies 
-z / .oyu4oj 


9A eo 
zo.oy 


9/1 ^7 
Z4. J / 


9/1 98 
Z4.Z5 


9/117 
Z4. 1 / 


s 

j 


offcenter clump on disk? 


SI 1 A87A7 
JJ. 105 /4 / 


97 AQOAQO 

-z / .oyu4yu 


^ 90 ss 
j> zy.oj 


9A 1 

zo. iy 


9S AS 
ZJ.D5 


9S S9 
Zj. JZ 


9 
z 


double clump 


SI 1 S3A 1 8 

JJ. IJJtlO 


97 AQ 1 31 A 

-z / .oy i j 10 


>, in no 


9A 70 
ZO. Iy 


9A A9 


9A 1 3 
ZO. 1 J 


A 
4 


several fuzzy 


SI 1A9898 
JJ. 14ZOZO 


97 AQ 1 A0Q 

-z / .oy iouy 


90 S3 

zy.oj 


97 1 
Z. 1 .iy 


9A QS 

zo.ys 


97 09 
Z / .UZ 


A 
4 


several fuzzy 


SI 1773A3 
jj.l / / joj 


97 A0 1 A77 

-z /.oyio / / 


9Q 7s 

> zy. / j 


9 A 70 
Z0. IV) 


9S 7A 
Zj. /4 


9S 3A 
Zj. j0 


1 


clump edge on 


S3 1A33SA 

JJ. It-JJJ't 


97 AC) 17 AO 
-z / .oy i /oy 


30 3A 


9A 87 
ZO.o / 


26 44 


9A 38 
ZO. JO 


| 


LlUIlip 


S3 1 98AS8 
Jj. 1Z5050 


97 AQ9 1 A3 

-z / .oyz 14 J 


98 on 
zo.yu 


9A 83 
ZO.oJ 


9A QA 

zo.yo 


9A AS 
ZO.OJ 


s 

J 


faint clumpy disk? 


S3 1 39 AAA 
JJ. IjZD^tt 


97 AQ9 1 07 

-z / .oyz iu / 


>, 90 3A 

> zy.jo 


97 1 
Z / . 1 V) 


97 1 S 

Z / . 1 j 


97 1 9 
Z / . 1Z 


1 
1 


faint clump 


S3 179Q3S 

jj.l /zyj j 


97 AQ9S1 9 
-z / .oyzj 1Z 


30 17 

J> jU. 1 / 


9A QS 

zo.yj 


9 A 30 
ZO. JU 


9A 08 
ZO.Uo 


9 

z 


faint double 


S3 1 S3738 
JJ. 1 JJ /jo 


97 AQ97Q7 

-z / .oyz /y / 


90 1A 

zy . 1*+ 


9A 7A 
ZO. /'f 


9S QS 

zj.ys 


9S 7Q 
Zj. / y 


1 

1 


faint clump edge on (tadpole?) 


S3 1A7AA7 
J J . 14 / 44 1 


97 AQ3S1 1 

-z / .oy j j 1 1 


^ 90 Q8 


9 A 80 
ZO.OU 


9S AQ 

zj.oy 


9S 38 
Zj. jo 


9 
Z 


double clump 


53.144209 


-27.693704 


29.86 


26.57 


26.24 


26.17 


1 


faint clump 


53.161386 


-27.694332 


27.15 


24.64 


24.31 


24.09 


2 


chain/three clumps 


53.134128 


-27.694236 


29.45 


27.29 


26.53 


26.14 


4 


fuzzy 


53.168314 


-27.694244 


28.74 


26.89 


26.55 


26.03 


4 


fuzzy 


53.174811 


-27.694860 


29.50 


25.51 


25.08 


24.84 


1 


clump (faint tail?) 


53.159129 


-27.695401 


> 30.07 


26.50 


26.09 


25.82 


4 


fuzzy 


53.141777 


-27.695409 


> 30.19 


26.76 


26.46 


26.39 


4 


faint fuzzy edge on 


53.176046 


-27.695456 


> 30.24 


27.22 


26.78 


26.56 


1 


faint clump 


53.161879 


-27.695566 


> 29.71 


25.89 


24.99 


24.76 


3 


tadpole 


53.156940 


-27.695521 


30.05 


27.19 


26.77 


26.76 


1 


faint clump close to face on spiral 


53.134625 


-27.695697 


> 29.11 


26.02 


25.65 


25.68 


4 


fuzzy 



TAB LE 2 — Continued 



aj2000 (degrees) 


<5j2O0O (degrees) 


U 


B 


V 


i 


Class 


Remarks 


53.130320 


-27.695894 


> 28.80 


25.50 


25.04 


24.81 


3 


tadpole 


53.166073 


-27.695896 


29.41 


25.96 


25.77 


25.76 


1 


fuzzy 


53.159225 


-27.695951 


29.04 


26.95 


26.91 


27.21 


1 


faint clump 


53.132461 


-27.696248 


> 29.16 


26.70 


25.76 


25.55 


3 


double clump 


53.143279 


-27.696576 


29.80 


26.99 


26.54 


26.51 


1 


fuzzy close to face on spiral 


53.138781 


-27.696842 


29.28 


26.90 


26.16 


25.60 


1 


fuzzy extended 


53.133936 


-27.696977 


28.51 


26.03 


25.96 


25.55 


1 


fuzzy extended 


53.151536 


-27.697612 


27.68 


25.60 


25.16 


24.87 


3 


tadpole 


53.136554 


-27.697853 


28.24 


27.08 


27.28 


27.36 


4 


fuzzy 


53.145437 


-27.697987 


26.50 


24.67 


24.46 


24.43 


3 


tadpole? 


53.173859 


-27.698226 


> 29.78 


26.40 


26.17 


25.65 


4 


faint fuzzy 


53.152117 


-27.698154 


28.77 


26.31 


25.97 


25.86 


3 


tadpole? 


53.174661 


-27.698140 


28.40 


26.58 


26.41 


26.39 


4 


fuzzy 


53.136690 


-27.698376 


28.87 


27.09 


27.04 


27.04 


4 


fuzzy 


53.160346 


-27.698494 


29.34 


26.98 


26.18 


26.05 


4 


fuzzy 


53.133898 


-27.698699 


27.34 


25.87 


25.71 


25.40 


4 


fuzzy 


53.149444 


-27.698816 


29.62 


27.29 


27.20 


26.96 


1 


faint clump 


53.182675 


-27.698849 


29.16 


26.41 


25.86 


25.59 


1 


clump 


53.162550 


-27.698872 


29.10 


26.68 


26.50 


26.72 


4 


fuzzy, fuzzy companions 


53.162513 


-27.699074 


29.20 


27.03 


26.64 


26.49 


4 


fuzzy, fuzzy companions 


53.178276 


-27.699784 


29.03 


26.95 


26.85 


26.81 


4 


fuzzy with fuzzy companion 


53.172033 


-27.684078 


28.57 


26.31 


25.78 


25.47 


4 


fuzzy 


53.180185 


-27.699950 


27.80 


25.90 


25.69 


25.56 


2 


faint double clump/tadpole? 


53.137583 


-27.700109 


29.01 


26.01 


25.46 


25.09 


1 


clump (spheroid) or faint star? 


53.169918 


-27.683871 


28.82 


26.61 


25.89 


25.60 


4 


fuzzy 


53.138746 


-27.700467 


> 29.76 


26.05 


25.26 


25.00 


2 


four clumps 


53.172890 


-27.683628 


> 30.20 


27.19 


26.77 


26.32 


1 


faint clump 


53.173737 


-27.700799 


29.27 


26.83 


26.67 


26.91 


4 


fuzzy 


53.172186 


-27.683993 


> 30.07 


27.25 


26.77 


26.30 


4 


faint fuzzy, close to clump 


53.143580 


-27.683502 


> 28.12 


25.72 


25.40 


25.24 


4 


fuzzy 


53.169913 


-27.669889 


> 30.22 


26.60 


26.17 


25.77 


1 


clump 


53.161684 


-27.682908 


28.35 


26.69 


26.70 


26.96 


1 


clump 


53.163174 


-27.669346 


27.54 


25.20 


24.80 


24.67 


5 


clumpy disk 


53.161053 


-27.672755 


29.26 


27.02 


26.41 


26.23 


1 


clump 


53.140722 


-27.682237 


> 28.08 


26.79 


26.86 


26.53 


4 


faint fuzzy 


53.166505 


-27.672539 


29.80 


27.19 


26.73 


26.54 


4 


fuzzy 


53.160749 


-27.671728 


27.15 


25.48 


25.23 


25.08 


4 


fuzzy 


53.164593 


-27.674991 


> 29.94 


26.48 


25.94 


25.89 


4 


fuzzy 


53.158301 


-27.680860 


> 29.15 


26.50 


25.70 


25.33 


4 


fuzzy 


53.169616 


-27.672793 


28.89 


26.71 


26.06 


25.62 


3 


tadpole 


53.158060 


-27.673707 


28.22 


25.78 


25.42 


25.19 


2 


three clumps, chain 


53.161286 


-27.679136 


28.91 


26.94 


27.00 


27.11 


1 


clump 


53.158157 


-27.673807 


27.09 


25.05 


24.83 


24.69 


2 


double clump 


53.162497 


-27.671140 


28.45 


26.32 


25.95 


25.85 


5 


fuzzy edge on disk 


53.165407 


-27.674053 


30.12 


26.92 


26.76 


26.57 


4 


faint fuzzy 


53.169762 


-27.674803 


28.17 


27.06 


27.28 


27.27 


4 


fuzzy 


53.152429 


-27.677562 


27.89 


26.39 


26.25 


26.16 


4 


fuzzy 


53.156011 


-27.676234 


28.80 


25.11 


24.83 


24.80 


3 


tadpole? or peculiar disk 


53.157810 


-27.676274 


28.84 


26.82 


26.62 


26.43 


4 


fuzzy 



TABLE 3 

Morphological classification of U band dropouts made using visual 

EXAMINATION OF GOODS/ACS Z-BAND IMAGES. THE RELATIVE FRACTIONS OF 
EACH CLASS ARE LISTED. 



Morphological class Fraction 



Single Clump 31% 

Double or multiple clump 13% 

Clump with tail 7% 

Extended and fuzzy 38% 



Fuzzy, extended with disk-like structure 11% 



